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Nomenclature 

a = absorption coefficient of semitransparent coating, 

„-i 

m 

c 0 = speed of electromagnetic propagation in vacuum, 

m/s 

G = flux quantity 2 (q, + ql), W/m 2 ; G L = G/oT 4 g[ in 
large-frequency spectral band 

h t , h 2 = convective heat transfer coefficients at boundaries, 
W/m 2 K; H = /i/oT 2 , 

K — extinction coefficient of semitransparent material, 
a + <r t , m~' 

k r , k,„ = thermal conductivity (coating, metal), W/mK; 

N r = k c /o-T 3 xl S c 

n = refractive index of semitransparent material 

q = heat flux, W/m 2 ; q r = radiative flux, W/m 2 ; 

q = qlcrTgi 

q q~ = isotropic radiative fluxes in and — x directions. 
Fig. 1, W/m 2 

T — absolute temperature, K; t — TlT xi 

T g u T g2 — gas temperatures on hot and cold sides of 
composite. Fig. 1, K 

T si , T s2 = temperatures of surroundings on two sides of 
composite. Fig. 1, K 

x = coordinate in composite, m: X = x/S c 

5 C , 8,„ = thicknesses of coating and metal, m 

s = emissivity 

Kcl = optical thickness of coating in large-frequency 

band, K L 8 C 

A, v = wavelength and frequency of radiation; 
v co — cutoff frequency 

p = diffuse reflectivity of interface from Fresnel 

equations 

cr = Stefan— Boltzmann constant, W/m 2 K 4 

ov = scattering coefficient in semitransparent coating, 

m“ ! ; Cl = o-JK 
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Subscripts 
c = coating 

i — inside interface at the hot side of coating. Fig. 1 

L — in spectral band where v > v co 

m = metal 

me = cold side of the metal wall 

mh = hot side of the metal wall 

o = outside interface at the hot side of coating. Fig. 1 

so = soot 

v - frequency dependent 

Introduction 

U SING thermal barrier coatings on combustor liners, tur- 
bine vanes, and rotating blades is important for reducing 
metal temperatures in current and advanced aircraft engines. 
Zirconia is a common coating material, and it is partially trans- 
parent to thermal radiation. Radiation becomes more signifi- 
cant as temperatures are raised for higher efficiency in ad- 
vanced engines. Calculations are often made with radiation 
effects neglected inside the coating. The effect of radiation is 
illustrated here, where an analytical procedure is provided by 
using the two-flux method for the radiative contribution. 

In Ref. 1 a detailed study was made of ceramic thermal 
barrier coatings for diesel engines. In Ref. 2 a two-flux analysis 
was developed for radiation in semitransparent multilayer 
composites. These references provide the basis for the present 
analysis where illustrative solutions are obtained for typical 
conditions in an aircraft engine. The formulation and solution 
of the exact spectral radiative transfer equations including 
large scattering, as is characteristic of zirconia, are rather com- 
plicated. For plane-layer composites, the two-flux equations 
were shown in Ref. 2 to give accurate results for gray and 
spectral layers; hence, the two-flux method is used here to 
provide a simplified method. 

Analysis 

A composite is first considered of a thermal barrier coating 
on a metal wall with the external surface of the coating cov- 
ered with a thin opaque layer of soot (upper part of Fig. 1). 
For an advanced engine, the combustion chamber pressure is 



Fig. 1 Geometry and nomenclature for a thermal barrier coating 
on a metal wall with and without soot on the exposed surface of 
the coating. 
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high enough that the combined gas and soot radiation can be 
assumed to provide a black environment at T sl . At the outer 
surface of the soot layer, the total heat flux consists of con- 
vection and radiation. 


at T s2 , such that radiative loss is by transfer between parallel 
plates. For e K , large surroundings at T s2 are assumed: 

<?,o. = h 2 {T mc - T s2 ) + On or e„)cr(n= - T , 2 ) (7a) 


<?,o, = h y {T gX - T so ) + e so ct(T 4 s1 - K) (1) 


e|| — ^mc), tmc 


(7b) 


The soot layer is very thin, and so temperature variation 
through its thickness is neglected. Within the semitransparent 
thermal barrier coating, the total heat flux is the sum of con- 
duction and radiation, 

, d7’ t (x) r 

<7,0, = ~K — + J q„ r {pc) dv ( 2 ) 

The q vr is eliminated by using the two-flux equation 3 (1/ 
K„)(dG„/dx) dv = —3 qj(x) dv, and the result is integrated in 
the coating from x — 0 to x, to yield the total heat flux as 


Two-Band Form of the Radiative Relations 

For a zirconia thermal barrier coating, the semitransparent 
region extends to a cutoff wavelength of about 5 p.m, beyond 
which the extinction coefficient becomes quite large. 1 ' 5 From 
the spectral property variations of zirconia, a two-band ap- 
proximation is reasonable, where the band for large wave- 
lengths (small frequencies) is considered opaque. 

To obtain the temperature distribution in the thermal barrier 
coating, Eq. (3) is integrated over the semitransparent band 
(large frequencies) and the result rearranged to give in dimen- 
sionless form 


q tot x = k c {T c (0) ~ T c (x)] + 


;L 


OX 0) ~ GJjx) 
K v 


dv (3) 


tJX) = / c (0) - — 


q**X 


cm - g l (X) 

3k cL 


( 8 ) 


The two-flux equation for G„ is a second-order differential 
equation that contains the local blackbody emission, 3 n 2 e vb [v, 
T(x)} 

d dv = (1 - £l v )[4n 2 e„ b (x) dv - G p (x) dv] (4) 

3 K v dx 

where 4 e^iv, T) dv = 27tC,v 3 dv/cl [exp(C 2 v/c 0 7’) — 1], ( C x 
and C 2 are the blackbody radiation constants). Two boundary 
conditions are required for solving Eq. (4). At the interface (x 
= 0) of the soot and the semitransparent coating, the radiative 
flux within the coating in the positive direction is composed 
of emission and reflection by the soot, qt r (Q) = e so n 2 e^(0) + 
(1 — s so )<?^(0). The two-flux relations, 3 qt r = (l/2)[(G„/2) + 
< 7 „ r ] and q~ r = (l/2)[(G„/2) — q vr ] are used to eliminate q* r and 
qZr, and the two-flux equation following Eq. (2) is used to 
eliminate q ur in terms of dG„/dx. This yields the boundary 
condition for G u at x — 0: 


dGp 

dx 


3 


2(2 - e so ) 


[G„(0) - 4n 2 <? l *(0)] 


(5a) 


A similar derivation gives the boundary condition for G„ at 
x = S,: 


To obtain G L (X), Eq. (4) is solved in the form 

= -3k c \(1 - CL L ){4n 2 tXX) A F[tXX)] - G L (X)} (9) 
dX 

where F(t) is the fraction of blackbody energy in the frequency 
range v = v co — °° for emission at t; F(t) was evaluated from 
the series form in Ref. 4, Chap. 2. For Eq. (9), there are two 
cases for boundary conditions: 1) an opaque soot coating on 
the zirconia, and 2) the zirconia surface is clean. For an opaque 
soot coating, Eq. (5a) applies at X = 0 

= - {G L {. 0) - 4rt 2 r(0) 4 F[/(0)]} (10a) 

x=o 2(2 £>so) 

When the zirconia surface is clean, q to , consists of radiation in 
the semitransparent spectral band combined with convection 
and radiative exchange at the surface in the opaque band. Us- 
ing the two-flux equation following Eq. (2) for the radiative 
transfer in the semitransparent band, gives 

= 3/c<J//,[l - / c (0)] + (1 - p o )0Ul - Ffc.)] 

x=o 

- tx 0) 4 { 1 - F[r c (0)]}) - q m ] (10b) 


d G l 
dX 


d Gl 
dX 


d G v 


- 4 " 2e *^>] ( 5b > 

2(2 - e mh ) 


If there is no soot coating, so that the thermal barrier coating 
has a clean external surface (lower part of Fig. 1), the bound- 
ary condition for G„ at the external boundary of the coating 


For a clean thermal barrier coating, Eq. (5c) gives 


GX 0) = 4 



rf.F(/„) 


2 1 + p, d G L 


3k cL 1 - p, dX 




(10c) 


At the interface between the coating and the metal wall, Eq. 
(5b) applies in a similar form to Eq. (10a). 


G„(0) = 4 


1 ~ Po 

1 ~ pi 


e vb (v, T sl ) 


+ 


2 1 4- pj d G, 

3 K v 1 — pi dx 


-V — 0 


(5c) 


At the interface between the semitransparent coating and 
metal, continuity of temperature gives T C (5 C ) = T m (8J = T mh . 
For the opaque metal wall, heat is transferred only by con- 
duction: 


<7,o, = k m [(T mh - r mc )/S„,] (6) 

At the cooled side of the metal wall, heat flow is by con- 
vection to the cooling gas, and by radiation to the surrounding 
structure. Since radiative transfer is appreciable, two limiting 
cases can be considered. For % there is an adjacent metal wall 


Solution Method 

An iterative solution was used. Since t(X) is unknown, the 
opaque heat conduction solution was used as a first guess for 
t c (X) and < 7 ,o t . When the coating surface is clean, dG L (0)!dX 
was found from Eq. (10b), and G^O) from Eq. (10c). Equation 
(9) for GXX ) was solved using Runge-Kutta numerical inte- 
gration. The temperature distribution in the coating was then 
obtained from Eq. (8). By using the temperature f t .(l) = t mh , 
Eqs. (6) and (7) in dimensionless form were solved simulta- 
neously for t mc and a new value of q, ot . The dG t (l)/dX from 
Eq. (5b), as in the form of Eq. (10a), was evaluated and com- 
pared with dG L ( 1 )/d X from the Runge-Kutta solution. The 
new t c (X) and < 7 , ot were adjusted by using a small fraction of 
the change from the previous iteration, and the adjustments 
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included making dG L (\)ldX from the dimensionless form of 
Eq. (5b) converge toward the Runge-Kutta solution value. At 
convergence, agreement was obtained to at least three signif- 
icant figures. The Runge-Kutta solution was checked by using 
20, 40, and 80 increments; no change to at least four significant 
figures was obtained by using 80 increments. For typical zir- 
conia thermal barrier coatings the absorption optical thickness 
is small and the scattering albedo is close to one. For these 
conditions the shooting method used for solving the two-point 
boundary value problem for G L converged without difficulty. 

The iteration method changes somewhat when the zirconia 
coating has an opaque soot layer on its exposed surface. The 
previous results for clean opaque and semitransparent coatings 
were used to estimate trial t c (X) and Gz.(0) values. With f so = 
t c ( 0), < 7 IO , was obtained from Eq. (1) in dimensionless form, and 
dG L (0)/dX from Eq. (10a). The iteration then proceeds as de- 
scribed for the coating with a clean surface. 



Results and Discussion 

Illustrative temperatures are shown in Fig. 2 for a zirconia 
coating on a high alloy steel combustion liner. The solid lines 
are for an uncoated metal wall; the lower line is for oxidized 
metal ( e m = 0.6) on both sides, and for the upper line the metal 
is also covered with soot on the combustion side (e so = 0.99). 
In both instances the temperatures are above the metal melting 
point, and so a thermal barrier coating must be used. A limiting 
calculation is for a zirconia coating assumed to be opaque 
(dashed lines). For a 1 -mm- thick coating the metal temperature 
is substantially reduced, and having soot on the coating has 
only a small effect as the coating thermal conductivity is 
mostly controlling. If the coating can be kept clean, there is a 
further benefit with the coating semitransparent. The high scat- 
tering of zirconia reflects away much of the incident radiation 
and temperatures in the zirconia are considerably decreased 
(lower dot— dash line), with the hot side of the metal reduced 
about 20 K. If, however, the zirconia becomes covered with 
soot, the metal temperatures are increased. The incident radi- 
ation is absorbed by the soot on the zirconia surface, and the 
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Fig. 2 Combustor liner wall temperature distributions for oxi- 
dized metal without a coating, with an opaque thermal barrier 
coating, and with a semitransparent thermal barrier coating with 
and without soot on its exposed surface. On the cooled side, ra- 
diation is to large surroundings, s x in Eq. (7a). Parameters (units 
are in the Nomenclature): A, = 250, h 2 = 110, k c = 0.8, k„, = 33, 6 C 
= 1<T\ 6 m = 0.794 x 10"\ n = 1.58, a = 30, and ar s = 10 4 for A < 
5 jum, r sl = T gl = 2000, and T s2 = T g2 = 800. 


Position within composite, x, mm 

Fig. 3 Turbine blade wall temperature distributions for clean 
metal without a coating, with an opaque thermal barrier coating, 
and with a semitransparent thermal barrier coating with a clean 
exposed surface (no radiative loss on interior side). Parameters 
(units are in Nomenclature): h l = 3014, h 2 = 3768, k c — 0.8, k m = 
33, S c = 0.25 X 10~\ S,„ = 0.762 x 10 \ n = 1.58, a = 30, and 
tr s = 10 4 for A < 5 pm, T sl = T gl = 2000, and T s2 = T g2 = 1000. 

semitransparency of the zirconia acts like increased thermal 
conductivity; hence, an opaque coating is more effective. 

Illustrative results are shown in Fig. 3 for a 0.25-mm-thick 
zirconia coating on a turbine blade. The inside of the blade is 
assumed isothermal, and so internal radiative heat loss is omit- 
ted. Two cases are considered for the exterior: 

1) The blade is surrounded by other cooled blades so that 
there is no external radiative heating. 

2) The blade is exposed to blackbody surroundings at T s i 
from the combustor. 

The zirconia surface is assumed to be clean. Temperatures 
are shown for bare metal, an opaque zirconia coating, and a 
semitransparent zirconia coating. With or without external ra- 
diative heating, the relative positions of the T(x) are the same; 
semitransparency yields higher metal temperatures than an 
opaque coating. For the higher convective heat transfer coef- 
ficients on a turbine blade compared with a combustor liner, 
the radiative effects are relatively too small for scattering to 
have sufficient benefit by blocking incident radiation for hot 
surroundings. Semitransparency increases the effective coating 
thermal conductivity for heating primarily by convection, and 
metal temperatures are increased. These results illustrate how 
the analysis can be used to examine the effect of semitrans- 
parency for various parameters; semitransparency is found to 
be beneficial for limited conditions. 
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